Spin-orbit coupling and g-factor of X-valley in cubic GaN 
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We report our theoretically investigation on the spin-orbit coupling and g-factor of the X-valley 
in cubic GaN. We find that the spin-orbit coupling coefficient from sp 3 d 5 s* tight-binding model is 
0.029 eV-A, which is comparable with that in cubic GaAs. By employing the k • p theory, we find 
that the gr-factor in this case is only slightly different from the free electron g-factor. These results 
are expected to be important for the on-going study on spin dynamics far away from equilibrium in 
cubic GaN. 
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Due to the existence of the wide energy gap between 
the conduction band and the valence band, GaN has 
been proposed to be a promising candidate for many- 
electronic applications, such as the solid-state ultravi- 
olet optical sources and high-power electronic devices. 1 
Recently, the discovery of the room-temperature ferro- 
magnetism in GaN based materials 2,3 highlights its pos- 
sible application in future spintronic devices^— Another 
outstanding property of GaN for realizing the spintronic 
devices is the extremely long spin lifetime^ because of 
the relatively weak spin-orbit coupling (SOC) compared 
to the narrow-gap III-V compounds, such as GaAs and 
InAs£i2 

For a detailed understanding of the spin dynamics, the 
SOC is essential.— In cubic GaN, the inequality of the 
cation and anion in the crystal leads to the bulk inver- 
sion asymmetry, which results in the Dresselhaus SOC. 11 
Up to date, the investigation on spin properties in GaN 
is focused on the low energy case, where only the low- 
est valley, i.e., T- valley is relevant. Recently, Fu and Wu 
reported the Dresselhaus SOC coefficient of the T-valley, 
0.51 eV-A 3 , from the sp 3 s* tight-binding (TB) model. 12 
This result agrees with the later experiment by the time- 
resolved Kerr rotation measurement. 13 However, for the 
spin dynamics under the influence of high electric field 14 
or with spin pumping by high-energy laser^ electrons 
can be driven into the high valleys This multivalley 
correlation was proposed to be able to induce the charge 
GmrniS and spin Gunnel effects in GaAs. Zhang et al*£- 
investigated the spin dynamics under the high electric 
field in GaAs quantum wells and suggested that the spin 
Gunn effect can be hindered by the fast spin relaxation 
of upper valley (L-valley) and hot-electron effect. To our 
best knowledge, there has been no report on the multi- 
valley spin dynamics in GaN till now. Since the upper 
valley in GaN is the A-valley (no L-valley exists in cubic 
GaN), one may expect different multivalley spin proper- 
ties in this material. Therefore, the details of the SOC in 
GaN is required. In our previous work, the expression of 
the SOC of the A-valley in cubic III-V semiconductors 
was derived^ but the corresponding coefficient of GaN 
is still unavailable. In the present communication, we 



calculate this coefficient for further investigations on the 
spin dynamics in this material. Moreover, the g-factor is 
also required to take into account the effect of the exter- 
nal magnetic field. Therefore, we will also calculate the 
g-factor of the A-valley in GaN from the k • p theory. 

In order to obtain the splitting energy of SOC, one 
needs to calculate the band structure. One of the most 
widely used approaches for band structure calculation 
is the TB theory— From our previous work on GaAs, 
we find that the d-orbitals are important for the spin 
splitting of the high valleys i 19 ' 21 Therefore, we employ 
the sp 3 d 5 s* nearest-neighbor TB model with the SOC 
herei 2 - 2 - The parameters are taken from the work by Jancu 
et al.M 

The spin-orbit splitting of the conduction band around 
the X-valley is plotted as a function of the momentum 
along X —> K and X — > W directions in Fig.QJa). One 
finds that the splitting increases linearly with the mo- 
mentum in the small momentum regime with respect 
to the bottom of the X-valley. In the large momentum 
regime, this monotonic tendency can be violated. 

For the states close to the A-point, the splitting can 
be described by the effective SOC Hamiltonian, f2(k) • 
<t, with ft and <x denoting the corresponding effective 
magnetic field of the conduction band and the Pauli 
matricesi 19 ' 24 For the valleys lie in the [001]-direction, 
one obtains^ 



n(k) = /3(k x ,-k v ,0). 



(1) 



Here k represents the momentum measured from the bot- 
tom of the valley. Obviously, this term results in the spin- 
orbit splitting linearly depending on the momentum, i.e., 

AE = 2/3fe|| with k\\ = ^Jk% + ky being the magnitude of 

the transverse momentum. 

The SOC coefficient can be measured by /3(k) = 
AE/(2kn). The value of j3 is shown as a function of mo- 
mentum along X — > K direction in Fig.QJb). One can see 
that in the range of small momentum, (3 keeps a constant 
value 0.029 eV- A. However, when the momentum lies far 
away from the bottom of the X- valley, (3 decreases. In 
the X — ¥ W direction, f3 increases with increasing mo- 
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FIG. 1: (a) Spin-orbit splitting of cubic GaN around the X- 
valley along X — > K and X — > W directions, (b) The cor- 
responding SOC parameter /3 vs. momentum along X — s> K 
direction, (c) /3 vs. angle for the momentum lying in the x-y 
plane. Solid curve: A;=0.01; dashed curve: 0.1; dotted curve: 
0.2; and chain curve: 0.3 (27r/o), with a — 4.5 A, the lattice 
constant of cubic GaN.— 



mentum (not shown) as expected from Fig.[TJa). This 
phenomenon is due to the higher order corrections of the 
SOC terms. We now turn to the effect of the direction 
of the transverse momentum on the SOC coefficient. In 
Fig-He), w e show the anisotropic behavior of SOC coef- 
ficient, where the momentum lies in the x-y plane with 
8 being the angle between the momentum and ir-axis. 
One can see the angle dependence becomes remarkable 
for large momentum, which is also from the higher order 
correction of the SOC. Interestingly, one notices that the 
SOC coefficient of X-valley obtained here is comparable 
with that of GaAs, 0.059 eV-A,i£ which is very different 
from the situation of T- valley. For the T-valley, the Dres- 
selhaus SOC coefficient in GaN is much smaller than that 
in GaAs. 

We should point out that the d-orbitals are of criti- 
cal importance in determining the SOC coefficient of the 
X-valley in GaN. Specifically, we obtain /3 = 0.002 eV- A 
from sp 3 s* TB model parameterized by O'Reilly et al.^ 
which is one order of magnitude smaller than that from 
sp 3 d 5 s* TB model. As a comparison, we also calcu- 
late the SOC coefficient of the T-valley from sp 3 d 5 s* TB 
model and obtain 7 = 0.235 eV- A 3 , which is close to that 



from sp 3 s* model 7 = 0.508 eV-A 3 (Ref.[U). 

Now, we turn to figure out the <?-factor of the X- valley 
based on the k • p approach by following the approach 
given in Ref.l26l. In our calculation, we first include the 
conduction band {X\ c ) and the valence bands (^3,, and 
X$ v ) and neglect the contribution of the other remote 
bands. 

Similar to the L-valley case^S one can write the lon- 
gitudinal effective mass mi and the transverse one mt 



,26 



and 



m 
mi 



1 = 



2 \(X lc \ Pz \X 3v ) 
m E Xla - Ex 3 , 



m , 2 \(X lc \ Px \X 5v )\ 2 



m t 1 m E Xx 



E 



(2) 



(3) 



respectively. 

The longitudinal and transverse g-factors are given 
by26,27 



.9|| - .90 



2 5{X lc \ Px \X 5v ){X 5v \p y \X 1 



m {E Xlc - Ex 3v )(E Xlc 
l)/(E Xlc - E X3 , 
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and 
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2 S'(X lc \p y \X 5v )(X 3v \p z \X lc ) 
m (E Xlc - E X3v )(E Xlc - E Xsv ) 

(m /m t - l)(m /m; - 1) 
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-S' 
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(5) 



where S = 2i(X^ v \h z \X^ v ) and 5' = 2i(X^ v \h x \X 3v ) are 
the matrix elements of the SOC>2Z, p describes the mo- 
mentum operator, go is the g-factor of the free electron. 

From our TB calculation, we obtain the energy levels 
of the X-valley, i.e., E Xlc = 4.58 eV, E Xsv = -2.74 eV, 
and E Xsv — —6.98 eV. We estimate the matrix element 
of the SOC from the band splittings of the X$ v band and 
take 8' = 8 = 0.02 eV. With mi = 0.5mo and m t — 0.3mo 
(Ref.H), we obtain = 1.996 and g ± = 1.997 by taking 
go = 2. We should point out that such a small differ- 
ence of the g-factor in X-valley from go is consistent with 
the previous results in silicon^ and other III-V group 
compounds^ 

Finally, we would like to discuss the contribution of 
the remote bands. One may notice that the second con- 
duction band (X 3c with E X3c = 8.21 eV) lies close to 
the Xi c bandj» Since its symmetry is the same as that 
of the X 3v band, it can also contribute to m; and g±. 
Similar to Eq. @ , one calculates the correction of the 
longitudinal effective mass due to the X 3c band and 
finds that the condition (Xi c \p z \X 3c ) / (Xi c \p z \X 3v ) < 
VWxTc - E x lc )/(E Xlc - Ex 3v ) ~ 0.56 to guarantee the 
real condition m; < mo- As a rough estimation, we take 
(X lc \p z \X 3c ) = 0.5(X lc \ Pz \X 3v ) and {X 5v | h x \X 3c ) = 
(X 5v \h x \X 3v ). Then, it is easy to calculate the fi- 
nal perpendicular g-factor g' ± = go + (g± — .90) (1 — 



3 



0.5ry)/Vl - 0.2577 with v = {E Xlc - E X3v )/(E Xac - 
Exic) — 3.18 and g± given by Eq. (|5|). One obtains 
g' ± = 2.004, which indicates that the relation g' ± w go 
is still preserving by considering the correction from the 
X^ c band. 

In summary, we have studied the SOC and <?-factor of 
the X-valley in cubic GaN. By taking into account the 
possible effect of the d-orbitals on high valleys, we per- 
form our calculation with an sp 3 d 5 s* nearest-neighbor 
TB model. The spin splitting of the conduction band in 
the X-valley and the corresponding SOC coefficient are 



calculated. We find that the SOC coefficient of the X- 
valley in GaN is 0.029 eV- A, which is comparable with 
that in GaAs. In addition, we calculate the g-factor and 
find that the value is very close to that of the free elec- 
tron. These results are useful for understanding the spin 
dynamics far away from the equilibrium. 
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